On the self-consistent general relativistic equilibrium equations of neutron stars 
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We address the existence of globally neutral neutron star configurations in contrast with the 
traditional ones constructed by imposing local neutrality. The equilibrium equations describing this 
system are the Einstein-Maxwell equations which must be solved self-consistently with the general 
relativistic Thomas-Fermi equation and /3-equilibrium condition. To illustrate the application of 
this novel approach we adopt the Baym, Bethe, and Pethick (1971) strong interaction model of the 
baryonic matter in the core and of the white-dwarf-like material of the crust. We illustrate the 
crucial role played by the boundary conditions satisfied by the leptonic component of the matter 
at the interface between the core and the crust. For every central density an entire new family 
of equilibrium configurations exists for selected values of the Fermi energy of the electrons at the 
surface of the core. Each such configuration fulfills global charge neutrality and is characterized by 
a non-trivial electrodynamical structure. The electric field extends over a thin shell of thickness 
~ h/{m e c) between the core and the crust and becomes largely overcritical in the limit of decreasing 
values of the crust mass. 
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One of the fundamental issues in physics and astrophysics is the creation of an electron-positron plasma in overcriti- 
cal electric fields larger than E c = m 2 c 3 /(eft) (see [l[ and references therein). Basic progress toward the understanding 
of the thermalization process of such a plasma have been achieved 0. The existence of such an electron-positron 
plasma has a central role in a variety of problems ranging from the acceleration process in gamma ray bursts (GRBs) 
[l[ to the sharp trigger process in supernova phenomena 0, ||J ■ This has motivated us to reconsider the standard 
treatment of neutron stars in order to find a theoretical explanation for the emergence of a wide variety of astrophys- 
ical situations involving such overcritical electric fields. In a classic article Baym, Bethe and Pethick |5[ presented the 
problem of matching to the crust in a neutron star a liquid core composed of N n neutrons, N p protons and N e electrons. 
After discussing various aspects of the problem they conclude: 'the details of this picture requires further elaboration; 
this is a situation for which the Thomas-Fermi method is useful.' In this letter we focus on relaxing the traditional 
condition of local charge neutrality n e = n p , which appears to have been assumed only for mathematical convenience 
without any physical justification. Instead, we adopt the more general condition of global charge neutrality N e — N p . 
The corresponding equilibrium equations then follow from self-consistent solution of the relativistic Thomas-Fermi 
equation, the Einstein-Maxwell equations and the /3-equilibrium condition, properly expressed in general relativity. 

The pressure and the density of the core are mainly due to the baryons while the pressure of the crust is mainly 
due to the electrons with the density due to the nuclei and possibly some free neutrons due to neutron drip (see e.g. 
Q). The boundary conditions determined by the matching of the electron distribution in the core with that of the 
electrons of the crust are fundamental for the self-consistent construction of the equilibrium configurations. 

We consider the case of a non-rotating neutron star with metric 

ds 2 = e v dt 2 - e x dr 2 - r 2 d8 2 - r 2 sin 2 9d<j) 2 , (1) 

where v and A are functions only of r. We assume units where G = ft = c = 1 and let a denote the fine structure 
constant. As usual we define the mass of the star M{r) by e" A = 1 — 2M/r + r 2 E 2 , and denote the Coulomb potential 
by V(r), which determines the electric field E = e~( v+x ^ 2 V' , where a prime indicates the radial derivative. The 
combined energy-momentum tensor of the matter and fields T^ v is given by 

Tp = diag(£ + £ em , -P - P em , -P + P e "\ -P + P em ) , (2) 

where £ cm = —P em = E 2 /(8tt), and £ and P are the energy density and pressure of matter. With all the above 
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definitions, the time-independent Einstein-Maxwell field equations read 

W = 47rr 2 £ + 47re A/ V eE(n p - n e ) , 
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In order to close the system of equilibrium equations, the condition of local charge neutrality n, 
traditionally imposed for mathematical simplicity. In this case the problem is reduced to solving only the Einstein 
equations for a Schwarzschild metric. When this condition is relaxed, imposing only global charge neutrality N e — N p , 
we need to satisfy the Einstein-Maxwell equations (J3])-©. In order to impose global charge neutrality as well as 
quantum statistics on the leptonic component, the general relativistic Thomas- Fermi equation must also be satisfied. 

The general relativistic electron Fermi energy is given by 
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where \i e = y/{Pi) 2 +m 2 e and Pf = (37r 2 n e ) 1 / 3 are respectively the chemical potential and Fermi momentum of 
degenerate electrons. From Eqs. ^ and |(7J) we obtain the general relativistic Thomas-Fermi equation 
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The /3-equilibrium condition is expressed by 
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In order to take into account the effect of the compression of the crust on the leptonic component of the core we 
solve the equilibrium conditions for the core within a Wigner-Seitz cell Q- The radius Rws of this cell determines the 
Fermi energy of the electrons of the core which has to be matched with the Fermi energy of the leptonic component 
of the crust. Global charge neutrality is specified by 
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From Eqs. © and (|10[) we can determine self-consistently the proton, neutron, electron fractions inside the core as 
well as the radius Rws of the Wigner-Seitz cell of the core 0]. 

The coupled system of equations consisting of the Einstein-Maxwell equations ((3J)— (jSJ) , the general relativistic 
Thomas- Fermi equation l[5]). the /3-equilibrium condition © along with the constraint (|10p needs, in order to be 
closed, an equation of state (EOS) for the baryonic component in the core and for the leptonic component of the 
crust. 

In order to illustrate the application of this approach we adopt, as an example, the Baym, Bethe, and Pethick 
(BBP) [| strong interaction model for the baryonic matter in the core as well as for the white-dwarf-like material of 
the crust. The general conclusions we reach will in fact be independent of the details of this model. 

At the neutron star radius r = R, all the electrodynamical quantities must be zero as a consequence of the global 
charge neutrality condition. Consequently, we have a matching condition with the Schwarzschild spacetime which 
imposes the boundary condition 
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The boundary conditions at the center correspond to M(0) = and the regularity condition to n e (0) = n p (0). 
From the /3-equilibrium condition ([5]), we can evaluate the central chemical potentials /^ e (0), /x p (0), and /z„(0), or 
equivalently, the central number densities n e (Q), n p (0), and n„(0) @. From Eq. ([7]) we also have the relation 
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Having determined the boundary conditions at infinity and at the center, we turn now to the matching conditions 
at the surface of the core. Following BBP [5(, the neutron profile at the core-crust interface is given by 

n n (z) = <™< + «°'' e - n c n rust )f(z/b) . (13) 

We have defined n c ° re = n n (R c ) and n^ rMSt = n n {Rws)- Here R c is the radius of the core defined as the point 
where the rest-mass density reaches the nuclear saturation density, i.e., p{R c ) — po ~ 2.7 x 10 14 g cm -3 5]. The 
function f(z/b) satisfies /(— oo) = 1, /(oo) = 0, where b ~ {n c ° re — n'£ ust )~ 1 ^ ~ l/m w As proposed by BBP, an 
appropriate choice for the function f(z/b) is the Woods-Saxon profile f(z/b) — (l + e z / b )~ 1 . The z-coordinate lines are 
perpendicular to the sharp surface separating two semi-infinite regions (core and crust) in the planar approximation 
Q; the neutron density approaches n c ° re as z — > — oo and n^ ust as z — » oo. 

The matching between the core and the crust occurs at the radius RwSi where we have V'(Rws) = by virtue of the 
global neutrality condition given by Eq. (fT0|) . and we also choose the value of the Coulomb potential V(Rws) = 0. 
From the electron chemical potential p e (Rws) at the edge of the crust, we calculate the corresponding neutron 
chemical potential fJt n (Rws) according to the BBP treatment. If p n (Rws) — m n > 0, neutron drip occurs. In this 
case, the pressure is due to the neutrons as well as to the leptonic component, so we have the inner crust (see Table 
H] and 0] for details). For larger values of the radii, i.e., for r > Rws the condition p n (r) — m n < is reached 
at pdrip — 4.3 x 10 11 g cm -3 and there the outer crust starts, with the pressure only determined by the leptonic 
component. If p n {Rws) ~ m n < 0, only the outer crust exists. 

For a fixed central rest-mass density p(0) ~ 9.8 x 10 14 g cm -3 and selected values of we have integrated the 
system of equations composed by the general relativistic Thomas-Fermi equation ([B"|). the /3-equilibrium condition ([5]), 
the Einstein-Maxwell equations ([3])- (0, with the constraint of overall neutrality (fTU|) . 

We found that although the electrodynamical properties of the core are very sensitive to the Fermi energy of the 
electrons (see Table H] for details), the bulk properties of the core like its mass and radius are not sensitive to the value 
of Eg . This is perfectly in line with the results of Ruffini ct al. in Q- 
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TABLE I: Results of the numerical integration of the BBP model for selected values of for p(0) ~ 9.8 x 10 14 g cm -3 . 
We show the mass and radius of the core M(R C ) and R c , the Coulomb potential at the center and at the core surface eV(0) 
and eV(Rc), the peak of the electric field in the core-crust interface E max , the rest-mass density at the edge of the crust 
p crus t = p(Rws), the mass of the crust M crust , and the inner and outer crust thickness A r c and A° c . 

Particularly interesting are the electrodynamical structure and the distribution of neutrons, protons, and electrons 
as the surface of the core is approached (see Fig. [1]) . It is interesting to compare and contrast these results with the 
preliminary ones obtained in the simplified model of massive nuclear density cores ■ The values of the electric field 
are quite close and are not affected by the constant proton density distribution assumed there. In the present case, 
the proton distribution is far from constant and increases outward as the core surface is approached. 

In conclusion, for any given value of the central density an entire new family of equilibrium configurations exists. 
Each configuration is characterized by a strong electric field at the core-crust interface. Such an electric field extends 
over a thin shell of thickness ~ 1 /m e and becomes largely overcritical in the limit of decreasing values of the crust 
mass and size (see Table U and Fig. [1}. 

These configurations endowed with overcritical electric fields are indeed stable against the quantum instability of 
pair creation because of the Pauli blocking of the degenerate electrons [l| . It is expected that during the gravitational 
collapse phases leading to the formation of a neutron star, a large emission of electron-positron pairs will occur prior 
to reaching a stable ground state configuration. Similarly during the merging of two neutron stars or a neutron star 
and a white-dwarf leading to the formation of a black hole, an effective dyadotorus Q will be formed leading to very 
strong creation of an electron-positron plasma. In both cases the basic mechanism which makes gravitational collapse 
depart from a pure gravitational phenomena is due to the electrodynamical process introduced in this letter. 

Finally, it is appropriate to recall that the existence of overcritical fields on macroscopic objects of M ~ M© and 
R ~ 10 km was first noted in the treatment of quark stars [§, [l(| HH, [13 • In that case the relativistic Thomas- Fermi 
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FIG. 1: Left column: the surface electric field in units of the critical field. Right column: the surface particle number density 
of neutrons (solid), protons (short-dashed), and electrons (long-dashed) normalized to the nuclear density for selected values 
of . First row: = 0.20m^, second row = 0.35m,r. 



equations were also considered. However, in all of these investigations, a hybrid combination of general and special 
relativistic treatments was adopted, resulting in an inconsistency in the boundary conditions (see [6j]). The treatment 
given here in this letter is the first self-consistent treatment of the general relativistic Thomas-Fermi equations, the 
beta equilibrium condition and the Einstein-Maxwell equations. Critical fields are indeed obtained on the surface of 
the neutron star core involving only neutrons, protons, and electrons, their fundamental interactions, and with no 
quarks present. 

While we were preparing our work an extremely interesting observational problematic has emerged from the Chandra 
observations of Cas A CCO [H E|. It is with a similar steadily emitting and non-pulsating neutron star that our 
theoretical predictions can be tested. In particular, the existence for each central density of a new family of neutron 
stars with a smaller crust than the one obtained when the local neutrality condition is adopted. 

Indeed, the existence of neutron stars with huge crusts, i.e., with both inner and outer crusts, is mainly a consequence 
of assuming no electrodynamical structure (i.e., assuming local neutrality) and of allowing electrons to have larger 
values of their Fermi energy Ef (see details in Q). It can also be demonstrated that no consistent solution of the 
Einstein-Maxwell equations satisfying the local n e — n p condition exists, even as a limiting case [fj. 
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